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ABSTRACT 

For the first time nonradial oscillations of superfluid nonrotating stars are self- 
consistently studied at finite stellar temperatures. We apply a realistic equation of 
state and realistic density dependent model of critical temperature of neutron and 
proton superfluidity. In particular, we discuss three-layer configurations of a star with 
no neutron superfiuidity at the centre and in the outer region of the core but with su- 
perfluid intermediate region. We show, that oscillation spectra contain a set of modes 
whose frequencies can be very sensitive to temperature variations. Fast temporal evo- 
lution of the pulsation spectrum in the course of neutron star cooling is also analysed. 
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1 INTRODUCTION 

Studying the pulsations of neutron stars (NSs) is very im- 
portant and actively developing area of research, since com- 
parison of pulsation theory with observations can potentially 
give val uable information abou t the properties of superdense 
matter (jAndersson et al.ll201lh . Yet, it is an extremely dif- 
ficult theoretical problem even for normal (nonsuperfluid) 
stars. Superfluidity of baryons additionally complicates the 
theory because in superfluid (SFL) matter one should deal 
with several independent velocity fields. As a consequence, 
the hydrodynamics describing pulsations of SFL NSs is much 
more complicated in_comparison to ordinary (normal) hy- 
drodynamics (e.g.. iGusakov 2007 ). 

In this letter we report on a substantial progress in 
modeling and understanding the nonradial oscillations of 
nonrotating SFL NSs in full general relativity. It should 
be noted that the nonradial oscillations of such stars 
have been intensively studi ed in the literature s tartin g 
from the seminal paper by iLindblom fc Mendelll ([199J). 
In particular, nonradial os cillations of general relativistic 
NSs were considered by IComer^ Langlois fc L in ( 199^ 



lAndersson. Comer, fc Langlois (120021') : lYoshida fc Led 
1 20031) : iLin. Andersson fc Comeij (|200a ). Because of the 
complexity of the problem, most of these papers used a 
simplified microphysics input, i.e., toy-model equations of 
state an d simplified models of baryon superfluidity (see, 
howe v er. iLin et al.l 120081: iHaskell . Andersson fc Passamontil 



I2OO9I : iHaskefl fc Anderssonll201oi ). Moreover, in all these 
studies SFL matter was treated as being at zero temper- 
ature, an assumption that is unjustified at not too low 
stellar temperatures and may lead to a quantitatively 



incorrect oscillation spe ctra (JGusakov fc AnderssonI I2OO6I : 
iKantor fc Gusakovi201ll : this letter). 

Here we improve on this by considering nonradial os- 
cillations of SFL neut ron stars at finite temper atures. We 
follow an approach of iGusakov fc Kantorl (I2OIII ) (hereafter 
GKll) which allows us to analyse oscillations of general rel- 
ativistic NSs employing realistic equation of state, density 
dependent profiles of nucleon critical temperatures and fully 
relativistic finite-temperature S FL hydrodynamics. 

As it was first found by iLindblom fc Mendelll (|l994 ). 
oscillation spectrum of a SFL NS consists of two distinct 
classes of modes, the so called normal and superfiuid 
modes. The frequencies of normal modes almost coin- 
cide with the oscillation frequencies of a normal star 
and hence are independent of temperature. The spec- 
trum of these modes is therefore very well studied 



in the lit erature (see, e.g., Thome fc Campolattard 



Ha ^ J 

Benhar. Ferrari fc Gualtieril 



iMcDermott. Van Horn fc HansenI 



2004 ). 
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19881 : 
the contrary, 

SFL modes can be very temperature-dependent. Using the 
approach of GKll they can be decoupled from the normal 
modes and studied separately. Since rad ial SFL modes have 
alread y been thoroughly analysed in iKantor fc Gusakovl 
(|201l!), here we focus on the nonradial SFL modes. In what 
follows, the speed of light c = 1. 



2 BASIC EQUATIONS 

In this section we briefiy discuss the equation describing 
SFL oscillation modes [Eq. ([S])]. The detailed derivation of 
this equation can be found in GKll. We consider a NS with 
nucleonic core and assume that both neutrons (n) and pro- 
tons (p) can be superfluid. Following GKll, we introduce 

where 



the baryon current density j'f-^s 



-J(n)+J 



(P)' 
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j(€, = n^ "" + ^'fcii-rfe) 



(1) 



is the current density for particles i = n or p (e.g., 
iKantor fc Gusakovll201lh . Here and below the summation 
is assumed over the repeated nucleon index A; = n, p. In 
Eq. ([l} rii is the number density; u'^ is the four-velocity 
of "normal" liquid component; Wn-N is the four-vector that 
characterizes motion of superfluid neutron (k = n) or proton 
(k — p) component with respect to normal matter. Finally, 
the symmetric temperature-dependent matrix Yik (= Yki) is 
a relativistic analogue of the SFL entrainment matrix. Since 
electrons (e) are normal, their current density is jf,, = ricU^, 
where rio is the electron number density. The quasineutrality 
implies no — rip (for simplicity, we ignore possible admixture 
of muons in the core). 

In this letter we study small-amplitude (linear) oscil- 
lations of a nonrotating star being initially in hydrostatic 
equilibrium. Hence, for the unperturbed star one has u'^ — 
(e~''",0,0, 0) and w'f^^ = luj^ , = 0, while the metric is 
dS^ = —e^dt^ + e^dr^ + r^dil"^, where r and t axe the ra- 
dial and time coordinates, respectively; ^{r) and A(r) are 
the metric functions; and fl is the solid angle in a spherical 
frame with the origin at the stellar centre. For definiteness, 
we assume that the unperturbed matter in the stellar core 
was in beta-equilibrium, S^ = ^n — /ip — ^c = 0, where 
/id is the chemical potential for particles d = n, p, e. We 
also restrict ourselves to oscillations with vanishing electri- 
cal current, j'f ^ — jT, , = 0. The latter condition couples the 
SFL degrees of freedom. 



(p) 



-(rpn/rpp)< 



(n)- 



(2) 



As it was shown in GKll the interaction between 
the SFL and normal oscillation modes is controlled by 
the coupling parameter s, which is given by s = 
{nedP/dnc)/{nh dP/drih) where P(nb,no) is the pressure 
and Tib = Un + Up is the baryon number density. This pa- 
rameter is small for a wide set of realistic equations of state, 
\s\ < 0.05 (see fig. 1 in GKll), so that the approximation 
of completely decoupled SFL and normal modes (s = 0) is 
already sufficient to calculate the pulsation spectrum within 
an accuracy of a few per cent. In the s = approximation 
the quantities jt^-^., P, and the metric g^i, remain unper- 
turbed for SFL oscillation modes. This opens up a possi- 
bility to formulate an equation describing SFL modes that 
depends only on SFL degrees of freedom, i.e., on 10^,^ This 
equation follows from the energy-momentum conservation 
and potentiality condition for motion of SFL neutrons. For 
a nonrotating star it takes the form (see GKll) 

ILU {tJ,„YnkW(k)j - rihW^n)]) = ^c dj {S fJ.°° ) , (3) 

where 5fj,°° = Sfie'^'^ and the disbalance 5fi{n-b, ric) equals 



Sfi = 



1^/2 



23ne 



Ynk dn^ ^ 
ric rib dx>^ (*' 



Ynk 

1 

Wb 



"(k) 



(4) 



In Eqs. ([3]) and (0) j = 1, 2, or 3 is the space index, dj = 
d/{dx-'), 23 = d8^{n\^,ng,)/dne, and all perturbations are 
assumed to be oc exp(ia;f). Because Eq. ([3]) is linear, one 
can generally present (5/t°° as: 5/i°°(r, f2) = (5/ii(r) Yim(fi), 
where Yimiyi) is the spherical harmonic [notice that Sfii{r) 
does not depend on the index m, see Eq. ((5}]. Combining 



then Eqs. (I2|), Q, and Q one obtains the following equation 






ft® 



5/11. 



(5) 



Here prime means derivative with respect to r and h — 
e"/^ nil {fin Ub y), where y = rib Ypp/[fin (VnnVpp - "i^np)] - 1. 
Eq. ([5]) determines the eigenfrequencies of SFL oscil- 
lation modes and is valid in the region of the stellar core 
where neutrons are superfluid (hereafter, SFL-region). This 
equation depends on the internal stellar temperature T only 
through the parameter y = y{T /Tcn,T /Tcp) , where Tdir) 
is the profile of critical temperatures for particles i = n or 
p. The boundary conditions to Eq. ([5} are following. If neu- 
trons in the stellar centre are superfluid, the regularity of 
the solution requires S^i oc r' at r — >■ 0. If the outer bound- 
ary of the SFL-region coincides with the crust-core interface 
(where r = -Rcc), then the SFL current should not penetrate 
the crust. This condition implies Sfj,'i{Rcc) = 0. Finally, if 
T is so high, that the SFL-region does not spread all over 
the core, then its boundaries are determined by the condi- 
tion T — T'cn(r) (see Sec. |3]for more details). In that case, 
the regularity of the solution at such boundaries requires 



Sf^'i 



o^-'/2,„2 



UJ^5fll/{h'^). 



Thus, we reduce the problem of calculation of SFL 
modes to solving simple second-order differential Eq. ([5]). 
For any fixed multipolarity I the solution to Eq. ([Sjl consists 
of a set of eigenfrequencies cj;„ and eigenfunctions 5jj.in{r) 
which difi'er by the number of radial nodes n — 0, 1, 2, . . . 



3 MICROPHYSICS INPUT AND NS MODEL 

Prior to studying the SFL oscillations using Eq. ([Sjl, one 
has to specify equation of state (including profiles of Ten 
and Tcp) and construct a hydrostatic model of an un- 
perturbed star. In addition, one has to specify the pro- 
file of internal stellar temperature. High thermal conduc- 
tivity lead s to a rapid equilibration of T in t he NS core 
(see, e.g., iGnedin. Yakovlev fc PotekhinI I200H ). As a re- 
sult, the red-shifted internal temperature T°° = Te"^^ 
becomes almost const a nt. M oreover, as it was shown by 
iGusakov fc AnderssonI (|2006l ). for the SFL-region to be in 
hydrostatic and beta-equilibrium it must be in thermal equi- 
librium. Thus, in what follows we assume that r°° = const 
in the SFL-region. 

In the pr esent letter we employ the equation o f state 
suggested by lAkmal. Pandharipande fc Ravenhalll (|l998r ) 
(APR). The coupling parameter s for such an equation of 
state is small, \s\ ~ 0.02. The density profiles of Ten and Tcp 
that we use here are shown in the left panel of Fig. [l] They 
do n ot contradict to results of m icroscopic calculations (see, 
e.g., iLombardo fc Schulzell200ll ) and are similar to the nu- 
cleon pairing models used to explain o bservations of the coo l- 
ing NS in Gas A supernova remnant (|Shternin et al.ll201ll ). 
For definiteness, all calculations are performed for a star of 
the mass IAMq and circumferential radius R = 12.2 km. 

Since T°° , but not T, is constant in the SFL-region, it is 
convenient to introduce the red-shifted nucleon critical tem- 
peratures T^ = e'^ " Tck {k — n, p) to analyse how the size 
of SFL-region changes with T°°. The functions T^{r) and 
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P [gtm ^] 



Figure 1. (color online) Left panel: Nucleon critical temperatures 
Tck versus density p (A: = n, p). Right panel: Red-shifted critical 
temperatures T^ versus radial coordinate r. 



T^ (r) are shown on the right panel of Fig.[T] The red-shifted 
proton critical temperature is high, T^{r) ~ 2 x 10® K, so 
that superfluid protons occupy the entire core almost im- 
mediately after the NS birth. The function T^{r) has a 
maximum 1^ ^ 5.1 x 10® K at r = r™'''' « 0.757?. Near the 
stellar centre the density varies slowly with r which results 
in a weak dependence of T^ on the radial coordinate. As the 
star cools down to T°° < TJ^, the SFL-region is formed, ini- 
tially, as a narrow spherical layer. Upon subsequent cooling 
the layer becomes wider and, for example, at T°° = 4x 10* K 
it is shown by the hatched region in the figure. As the tem- 
perature decreases further, the SFL-region extends to the 
crust and, eventually, at T°° = T;^{0) ^ 2 x 10* K it pene- 
trates the stellar centre. 
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4 OSCILLATION SPECTRA AND MODES 

Figure [5] presents normalized eigenfrequencies ijJi„ (in units 
oi CO = c/R ~ 2.5 X lO** s~^) versus internal temperature 
T°° for SFL oscillations of multipolarity / = 0, 1, 2, and 3. 
For each I we plot a set of oscillation modes that differ by 
the number of radial nodes n = (solid lines), n = 1 (dots), 
n = 2 (long dashes), n = 3 (long-short dashes), and n = 4 
(dashes). One sees that the higher the n the larger the uin- 
By the hatches in the bottom panel of the figure we show 
the SFL-region; as expected, the size of this region depends 
on T°° . The grey-shaded area corresponds to the crust and a 
region in the core where all neutrons are unpaired. A similar 
shaded area on the four upper panels shows temperatures 
T°° ^ T^ for which all neutron matter in the core is normal. 
In the latter case there are no SFL modes in NS. 

Before further discussing spectra in Fig. 2 it is con- 
venient to describe briefly Fig. 3 that presents eigenfunc- 
tions S^i„ (r) normalized to unity in the maximum. The solid 
lines correspond to radial oscillation modes (1 = 0), dotted 
and dashed lines describe dipole (/ = 1) and quadrupole 
{I = 2) modes, respectively. Each column in the flgure con- 
tains four panels which are plotted for the following temper- 
atures (from bottom to the top): T°° = 10*, 2 x 10*, 3 x 10*, 
and 4 x 10* K. For any of these temperatures we have flve 
panels in a row, which correspond to (from left to right) 
n = 0, 1, 2, 3, and 4 radial nodes of (5/ij„(r). The stellar 
regions where neutrons are normal, are shaded in Fig. [3] 

At T°° < T^ the size of SFL-region rapidly increases 



Figure 2. (color online) Eigenfrequencies a;;„ versus T°° for mul- 
tipolarities I = 0, 1, 2, and 3. For each I a set of curves is plotted 
with n = 0, 1, 2, 3 or 4. At T°° ^ Tcn(O) Ri 2 x 10* K (see the 
vertical short-dashed line) superfluidity occupies the stellar cen- 
tre. The bottom panel shows variation of SFL-region with T°° . 
For more details see the text. 



as the star cools down, whereas the eigenfrequencies are al- 
most temperature-independent for T°° > 3 x 10* K (except 
for the modes with n = 0, see Fig. 2 and discussion below). 
This is so because of compensation of two opposite tenden- 
cies: (i) expansion of the SFL-region and (ii) increasing of 
the local speed o f SFL sound «°^ = e~'^^'^V- h JB with de- 
creasing T°° (see lGusakov fc Anderssonll2006l for details on 
v^^). The tendency (i) leads to decreasing while (ii) leads to 
growing of the eigenfrequencies Luin. 

The qualitative behaviour of cj;„ and 5^in{r) at 3 x 
10* K< r°° < 1"^ follows from simple arguments. The pa- 
rameter h is small there, so that it is sufficient to retain only 
the terms oc /i^Mn Eq. dS]). As a result, for T°° > 3 x 10* K 
SiJ.in{r) will be almost independent of I (see Fig. [3|. Fur- 
thermore, at T°° > 4 X 10 K the function h can be ap- 
proximated as h X A [(TS - r°°) - B{r - r^rf] , where 
A and B are some constants depending on a NS model. 
Using this approximation, one can analytically solve Eq. 
((SI and find that eigenfunctions 5jj,i„{r) are proportional 
to the Legendre polynomials P„ while the eigenfrequencies 
ujin cx \/n{n + 1) and are independent of T°° (see Figs. [2] 
and |3] for r°= > 4 x 10* K). For the modes with n = 
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Figure 3. (color online) Eigenfunctions Sfii„ versus r for T°° = 10*, 2x10*, 3x10*, and 4x10* K (5 panels in a row for each temperature). 
The columns of panels are for n = 0, 1, 2, 3, and 4 radial nodes. Solid, dotted, and dashed lines correspond to multipolarities Z = 0, 1, 
and 2, respectively. For more details see the text. 




Iff 



10* 

t [years] 



Iff 



-200 -100 100 200 300 400 

t — ti [years] 



10 12 

t — ti [years] 



Figure 4. (color online) Upper panels: The oscillation frequencies U2n of quadrupole modes with n = 0,. . .,4 versus stellar age t. 
Vertical dashed lines show the stellar age t\ at which neutron superfluidity first appears in the stellar centre. Bottom panels: Variation 
of SFL-region with t. For more details see the text and captions to Figs. [2] and [3] 



this estimate gives toio = 0. As follows from Fig. [2] in that 
case the eigenfrequencies are indeed small but nonzero. To 
improve the estimate for cJio one should take into account 
the term depending on l{l + l)/r^ in Eq. ((5]). Because at 
T°° > 4 x 10** K Sfj,io{r) is almost constant (see top- left 
panel in Fig. |3}, one can put 5/iJo ~ and S^'/q « 0. Us- 
ing then Eq. ([5}, one obtains the following estimate for cjiq: 
LOio oc VU^ + 1) «'7C„"" oc ^Z(Z-H)Vl-r-/'Igg. This 
formula indicates that ujio depends only on T°° /1^ (but not 
on the size of SFL-region) and vanishes at T°° — 1^ (see 
Fig- El . It remains to note that for / = Eq. (O has a static 
solution Sfioo = 5/i°° — const and ixjqo — 0, which describes 
a star in hydrostatic and diffusive eq uilibrium (but not nec- 
essarily in beta-equilibrium, see Gusakov fc Anderssonll2006l 
for more details). It should also be stressed that approximate 
formulas for eigenfrequencies found above is not a special 



feature of our microphysical model. Rather, it is inherent in 
all NSs for which the maximum of T^{r) lies between the 
centre and the crust-core boundary. 

At r°° slightly exceeding T;sIo) the size of SFL-region 
rapidly expands to the stellar centre as the star cools down. 
The reason for that is weak dependence of T^ on r at r < 
0.1 ii (see the right panel of Fig. [!)). Since T°° ^ T^ for 
r < 0.1 J?, v^^ is close to there. This results in a noticeable 
decrease of a;;„ with n ^ 1 near r;,^(0) (see Fig. [2)|. In 
particular, at T°° = 2 x 10* K the corresponding oscillation 
modes are all localized in the vicinity of the stellar centre, 
where w° is small (see Fig. [Sjl. During subsequent cooling 
[for T°° < rc^(O)] the entire core is already occupied by the 
neutron superfluidity and increasing of v^ , especially in the 
central regions of the star, leads to growing of aj;„. When 
T°° drops below ~ 4 x 10^ K the eigenfrequencies cj;„ and 
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eigenfunctions 5jj,i„{r) become almost independent of T°°, 
since h (and v^ ) approach their asymptotes at r°° — 0. 

Let us now discuss the temporal evolution of oscilla- 
tion spectra during the star cooling. To simulate cooling 
we apply a slightl y upda ted version of the code discussed 
by iGusakov et alj (|2005l V The onset of neutron superflu- 
idity is accompanied by an abrupt acceleration of cooling 
due to the Cooper pairing neutrino emission process (see 
IShternin et alJ I2OIII for recent observational evidences of 
such cooling). As is demonstrated in Fig. 2] this process 
leads to a rapid evolution of oscillation spectrum. Upper 
panels in Fig. Ushow the eigenfrequencies i'2n = aj2n/(27r) of 
quadrupole modes (n = 0, . . . , 4) versus NS age t. The lower 
panels illustrate change of the SFL-region with t. On the 
left upper panel the functions i'2n{t) are plotted for the time 
interval from 200 to 2 x 10^ years. After the onset of neu- 
tron superfluidity (at t ~ 330 years) the oscillation spectrum 
rapidly changes and approaches its zero-temperature asymp- 
tote only at f > 10'' years. At the star age ti ~ 480 years 
superfluidity penetrate the stellar centre [this corresponds 
to T°° = T;^{0)]. The period of time |t - fi| < ii is accom- 
panied by a very fast evolution of the spectrum. In more 
detail the evolution is shown in the central panels which 
span an interval of approximately 600 years of NS cooling. 
The time t in these panels is counted from ti. The right 
panels in Fig. |4] show the only four- year episode of NS life in 
the vicinity of ii. One sees that even for such a small period 
of time the frequencies of some oscillation modes increase 
by more than 10%, whereas the frequency of the mode with 
four radial nodes changes noticeably on a time-scale of a few 
months (!). 



5 SUMMARY AND OUTLOOK 

In this letter we study, for the first time, nonradial oscilla- 
tions of SFL NSs at finite temperatures. Use of an approach 
developed in GKll enable us to solve this problem in full 
general relativity, employing the realistic equation of state 
(APR) and realistic, density-dependent profiles of nucleon 
critical temperatures (Fig.[T]). It is shown that equations de- 
scribing SFL modes can generally be reduced to the simple 
second-order differential Eq. ([5]). The eigenfrequency spec- 
trum (Fig. ^ and eigenfunctions (Fig. ^ for this equation 
are carefully analysed. It is demonstrated that dependence of 
eigenfrequencies lji„ on T°° is determined by two competing 
effects: (i) decreasing of ujin with expanding SFL-region and 
(ii) growing of coin with increasing of the local speed of SFL 
sound v" . These results agree with the conclusions made by 
Kantor and Gusakov (2011) for a radially oscillating NS. 

In addition, we examine the evolution of oscillation 
spectrum in the course of NS cooling (Fig. [4]). We find that 
acceleration of NS cooling soon after the onset of triplet 
neutron superfluidity leads to a very fast modification of 
the spectrum - the eigenfrequencies can vary dramatically 
on a time-scale of months. 

In the present work we completely ignore various dissi- 
pative effects in pulsating SFL NSs. Meanwhile, our results 
provide a strong basis to study these effects using the per- 
turbative scheme suggested in GKll. In particular, to de- 
termine the damping time r of some normal or SFL mode 
one could proceed in the following steps: (i) find the vec- 



tors u'^ and w'f^^ assuming s = and neglecting dissipation 
(see GKll and Sec. 2); (ii) use them to calculate the dis- 
sipative terms entering the SFL hydrodynamics (and hence 
the rate of change E of the oscillation energy E, see, e.g., 
iKantor fc Gusakov|[2oTII ) : (iii) calculate r as r = -E/{2E). 
The other important problem concerns gravitational ra- 
diation from pulsating SFL NSs. The radiation from normal 
modes can be accurately calculated already in the s = 
limit (and will be the same as that for a nonsuperfluid NS). 
In turn, SFL modes are not coupled with the metric in the 
s = limit, so that to calculate gravitational radiation from 
them one must use the first-order perturbation theory in s. 
Interestingly, this problem is equivalent to finding gravita- 
tional radiation from a normal NS experiencing oscillations 
under an action of a small external force ex s. An intensity 
of such radiation will be reduced by a factor of s'^ ~ 10~^ 
in comparison to normal modes with the same E. Detailed 
studies of dissipation in SFL NSs are a very important task, 
e.g., for calculation of the instability windows of r- modes; 
we plan to address it in the near future. 
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